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Abstract 



This work is concerned with the formulation of a general framework for the analysis of 
meshfree approximation schemes and with the convergence analysis of the Local Maximum- 
Entropy (LME) scheme as a particular example. We provide conditions for the convergence in 
Sobolev spaces of schemes that are n-consistent, in the sense of exactly reproducing polynomials 
of degree less or equal to n > 1 , and whose basis functions are of rapid decay. The convergence 
of the LME in Wi^^{il.) follows as a direct application of the general theory. The analysis 
■^r ■ shows that the convergence order is linear in /i, a measure of the density of the point set. The 

^' I analysis also shows how to parameterize the LME scheme for optimal convergence. Because 

of the convex approximation property of LME, its behavior near the boundary is singular and 
requires additional analysis. For the particular case of polyhedral domains we show that, away 
Cu ' from a small singular part of the boundary, any Sobolev function can be approximated by means 

of the LME scheme. With the aid of a capacity argument, we further obtain approximation 
results with truncated LME basis functions in H^{n) and for spatial dimension d > 2. 

00 ■ 1 Introduction 

m ■ 

^^ ' Meshfree approximation schemes (cf., e. g., [13] for a review) are advantageous in a number of 

areas of application, e. g., those involving Lagrangian descriptions of unconstrained flows (cf., e. g., 
|16j for a representative example) where methods based on triangulation, such as the finite-element 
method, inevitably suffer from problems of mesh-entanglement. The present work is concerned with 
the formulation of a general framework for the analysis of meshfree approximation schemes (cf., 
e. g., ^7\ for representative past work) and with its application to the Local Maximum-Entropy 

K^ \ (LME) scheme as an example. By way of conceptual backdrop, we may specifically envision time- 

^ ■ independent problems for which the solutions of interest follow as the minimizers of a functional 

^ F : X ^ W, where X is a topological vector space. General conditions for the existence of solutions 

are provided by the Tonelli's theorem (e. g., [10]). In this framework, an approximation scheme is 
a sequence X^ of subspaces of X, typically of finite dimension, defining a corresponding sequence 
of Galerkin reductions of F, 

\ -|-oo, otherwise. 
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An approximation scheme is then said to be convergent if it has the fohowing density property: For 
every u € X, there exists a sequence Uk S X^ such that \iuik-^oo Uk = u. The connection between 
density of the approximation scheme and convergence is provided by the foUowing proposition [7]. 

Proposition 1. Let X be endowed with two metrizable topologies S and T, with T finer than S. 
Let F: X ^M. he coercive in {X^S) and continuous in {X,T). Let Xk be a dense sequence of sets 
in {X,T) and let F^ he the corresponding sequence of Galerkin reductions of F. Then the sequence 
Fk T -converges to the lower semicontinuous envelope of F and is equicoercive in {X,S). 

We recall that F-convergence is a powerful notion of variational convergence of functionals that, 
in particular, implies convergence of minimizers. Thus, if the sequence F^ is equicoercive, then 
the minimizers of F are accumulation points of minimizers of Ff., i. e., if F^^Uk) = inf F^ then 
the sequence Uk has a subsequence that converges to a minimizer of F. We also recall that the 
topology T is finer than S, i. e., any converging sequence for T converges for S. In applications, T 
is typically a metric or normed topology and S the corresponding weak topology. 

It thus follows that, within the general framework envisioned here, the analysis of convergence of 
approximation schemes reduces to ascertaining the density property. Towards this end, in Section [3] 
we begin by analyzing meshfree approximation schemes that are n-consistent, in the sense of exactly 
reproducing polynomials of degree less or equal to n > 1, and whose basis functions are of rapid 
decay. Specifically, for schemes subordinate to point sets possessing a certain geometrical regularity 
property that we term h-density, we prove a uniform error bound for consistent and rapidly-decaying 
approximation schemes. In addition, we show that the sets of functions spanned by consistent and 
rapidly-decaying approximation schemes are dense in Sobolev spaces. 

In Sections |4] and EJ we apply the general results of Section [3] to the Local Maximum-Entropy 
(LME) approximation scheme of Arroyo and Ortiz [2] (see also [3l [271 El H] ) • The LME scheme 
has been extensively assessed numerically over a broad range of test problems [21 HH [16], but a 
rigorous convergence analysis has been heretofore unavailable. The general theory of Section [3] 
readily establishes the density of the LME approximation spaces Xk in Wjq'^(O), cf. Section [H In 
particular, the analysis shows that the convergence order is linear in h, a measure of the density of 
the point set. These convergence rates and the corresponding error bounds are in agreement with 
the numerical results reported in [2] , and are comparable to those of the first-order finite element 
method (cf., e. g., [5]). Conveniently, the analysis also shows how to choose the LME temperature 
parameter so as to obtain optimal convergence. This optimal choice is in agreement with that 
determined in [21 [16] by means of numerical testing. 

The LME scheme is a convex approximation scheme in which the basis functions are constrained 
to take non-negative values. By virtue of this restriction, the LME scheme is defined for convex 
domains only. Consequently, its behavior near the boundary is somewhat singular and requires 
careful additional analysis. In Section[5l for the particular case of polyhedral domains we show that, 
away from a small singular part of the boundary, any Sobolev function can indeed be approximated 
by means of the LME scheme. Then, with the aid of a capacity argument we obtain approximation 
results with truncated LME basis functions in H^{Q,) and for spatial dimension d > 2. 

2 Prolegomena 

The open d-ball B{x,6) of radius 6 centered at x is the set {y € M*^: jy — x| < 6}. The closed 
d-ball B[x, 6) of radius 6 centered at x is the set {y G M : |y — 2;| < 6}. Given a set A C M , we 



denote by A its closure, and by dA its boundary. By a domain we shall specifically understand an 
open and bounded subset of M*^. Given a point set P € (M'^)^, we denote by conv(P) its closed 
convex hull [20], and by conv(P) the interior of conv(P). We recall that a d-sinaplex T C M'^ is the 
convex hull of d + 1 affinely independent points [20]. Given a bounded set ^ C M*^, its size Ha is 
the diameter of the smallest ball containing A. 

The following definitions formalize the notion of a point set P dVt that approximates a domain 
VL uniformly. 

Definition 1 (/i-covering). We say that a point set P C M. is an /i-covering of a set A C M , 
^ > 0, if for every x G A there exists a d-simplex Tx of size hx^ < h and with vertices in P such 
that X (z Tx- 

Definition 2 (/i-density). We say that a point set P C M"' has /i-density bounded by t > if for 
every x eW^, #{P n B{x, h)) < r. 

For a point set P C 0, with /i-density bounded by r, the following proposition bounds its number 
of points in rings of M . 

Proposition 2. Assume P C ^ has h-density bounded by t, for some h, t > 0. Then there is a 
constant c > that depends on t and d such that, 

#{Pr\ {B{x, th) \ B{x, {t - l)h))) < cf^-^, (2) 

Vx G il and integers t > 1. 

Proof. Let 

Ei = iyeR'^: dist {y, B{x, th) \ B{x, {t - l)/i)) < h\ 
= B{x, {t + l)/i) \ B{x, {t - 2)h) 



(3) 



and 



E2 = ly eM'^: dist {y, B{x, th) \ B{x, {t - l)h)) < 2h\ 
= B{x, (t + 2)h) \ B{x, {t - 3)h). 



(4) 



Then for every y G B{x, th) \ B{x, (t — l)h) there is a z ^ Z := Ei H hd 2Z such that \y — z\ < h 
and so 

B{x, th) \ B{x, {t - l)h) C y B{z, h). (5) 

On the other hand, z + [0, hd~2)^ c £'2 for all z e Z, and thus 



#Zh'^d~2 
Consequently, 



IJ z+[0,hd~^ 



2)'^ 



z^Z 



< \E2\ < \B{0,l)\h'' ({t + 2)" - {t - 3)") KchH"^'. (6) 



#{Pn{B{x,th)\B{x,{t-l)h))) <# (pn U Biz,h)] <cdiTt'^~\ 



(7) 
D 



3 Convergence Analysis of General Meshfree Approximation Schemes 

In this Section we analyze meshfree approximation schemes that are n-consistent and whose shape 
functions are of rapid decay. Specifrcahy, we prove a uniform error bound for consistent and 
rapidly-decaying approximation schemes. In addition, we show that the set of functions spanned 
by consistent and rapidly-decaying approximation schemes are dense in Sobolev spaces. 

Let il C M be a domain. By an approximation scheme {/, W, P} we mean a collection W = 
{wa, a G /} of shape functions and a point set P, both indexed by /. Given an approximation 
scheme {I,W,P}, we approximate functions u: f] — )■ M by functions in the span X of W oi the 
form 

Ul{x) = '^u{Xa)Wa{x), (8) 

provided that this operation is well defined. More generally, we shall consider sequences of approx- 
imation schemes {Ik,Wk,Pk} and let 

Uk{x) = ^ u{Xa)Wa{x), (9) 

ae/fc 

be the corresponding sequence of approximations to u in the sequence X^ of finite-dimensional 
spaces of functions spanned by W^. We note that, for simplicity, we assume that all functions are 
defined over a common domain Q. Depending on the approximation scheme, this assumption may 
implicitly restrict the type of domains that may be considered, e. g., polyhedral domains. The aim 
then is to ascertain conditions on the approximation scheme under which u^ — )• u in an appropriate 
Sobolev space iy'"'P(0). 

We recall the following definition of consistency of approximation schemes [23] . 

Definition 3 (Consistency). We say that an approximation scheme {I,W,P} is consistent of order 
n >0, or n-consistent, relative to a point set P if it exactly interpolates polynomials of degree less 
or equal to n within Q, i. e., if 

x- = ^x>„(x) (10) 

a&I 

for all multiindices a of degree \a\ < n. 

A simple binomial expansion shows that (jlOp can equivalently be replaced by 

^^Wa{x) = l, (11a) 

^Wa{x){xa - xY = 0, Vo G N'^, < |a| < n, (lib) 

in the definition of consistency. 

Consistency results in a number of identities involving the partial derivatives of the shape 
functions, which we record next for subsequent use (cf. [6]). 

Lemma 1. Let {I,W,P} he an approximation scheme. Suppose that W consists of C'^{Q) shape- 
functions that are nth-order consistent relative to P in Q. Let a, (3 he multiindices, with < |a| < n, 
< |/3| < r. Then, 

Y.D^u,^.)(x.-.r = \f 'f;'^- (12) 

— ; U, otherwise. 



Proof. We proceed by induction on |/3|. For /3 = 0, the identity (J12p follows directly from consis- 
tency. Let < m < r and assume eq. ()12p holds for all multiindices a, (5 such that < |/3| < m 
and < |q| < n. Let /3 and 7 be such that |/3| = tti and I7I = 1. Then, for a = we have, by 
consistency, 

Y, D^+^Wa{x) = D^+"' Y, ^a(x) = 0, (13) 

whereas for < \a\ < n we have, also by consistency, 

= D^Y. D^Wa{x){Xa -xr+Y DPwa{x){a ■ ^){Xa " x)""^ (^4) 

= {a-^)Y,DPwa{x){Xa-xr-^ 

Suppose that a 7^ /3 + 7. Then, from (fT2]) . 

YD^M^)ixa-xr-^ = 0. (15) 

Suppose, contrariwise, that a = /3 + 7. Then, also from (fT2]) . 



(a • 7) J^ D^«;a(x)(x, - xr-^ = (a . 7)(a - 7)! = «!, (16) 

whereupon p^ becomes 

^D''+^w;„(x)(xa-x)° = a!, (17) 

and ()12p holds for all multiindices /3 of degree m + 1. D 



We recall that the Taylor approximation of order r of a function u £ C^~^^(Q) at y G $7 is 

Tr{u){x,y) = Y ^D"^^i^){y - xT (18) 

and its remainder is 

Rr+i{x,y) = u{y)-Tr{u){x,y), (19) 

which turns out to be 

Rr+i{x,y)= V l-D''u{x + e{y-x)){y-xr, (20) 

\a\=r+l 

for some 9 € (0, 1). 

Functions in the span of a consistent shape- function basis satisfy the following multipoint Taylor 
formula (cf. [5l[6j). 



Proposition 3 (Multipoint Taylor formula). Let W he a C"'(0) shape- function set nth- order con- 
sistent relative to a point set P in $7, u G C^"'^(conv(J7)) and m = min{n,^}. Then, 

D'^Uiix) = D^U{X) + Y,Rm+liXa,x)D''Wa{x). (21) 

ael 

for all \a\ < min{m,r} and x E f]. 

Proof. The proof follows that of Theorem 1 of [6] . From the Taylor expansion of order m of u at 
X we have 

u{Xa)= Yl -o^D'^u{x){Xa-xf + Rn,+i{x,Xa). (22) 



|/3|<m 



whence it follows that 

D''ui{x) = Y,u{xa)D''wa{ 






J2\ Yl iD^uix){Xa-xf + R^+iix,Xa)\D''Wa{x) 
ael \|/3|<m / 

J2 i^M^) (5^^"w'a(x)(Xa-x)^ J +J2Rm+l{x,Xa)D' 
\B\<m \a£l ) aGl 



(23) 

Wa{x), 



and (|2ip follows from Lemma [TJ D 

We recall that a function / € C°°(M ) is said to be rapidly decreasing if [21] 



sup sup(l + |xn'''|(D"/)(x)| <oo (24) 

\a\<N xeW^ 



for all iV = 0, 1, 2, . . . , where jxp = ^ xt- 



"I- 
The next definition formalizes a polynomial-decay condition of the shape functions and their 

derivatives. 



Definition 4 (Approximation scheme with polynomial decay). We say that an approximation 
scheme {/, W, P} has a polynomial decay of order (r, s) for constants c > and h > if the basis 
W is in C^{Q), and 



sup sup sup 1 

\a\<rxen a£l \ 






h 



2^ 



h\'^\\D''wa{x)\ < c. (25) 



A sequence of approximation schemes {Ik,Wk,Pk} has a uniform polynomial decay of order {r,s) 
if there exists a constant c > and a sequence hj. —^ such that, for each k, {Ik,Wk,Pk} has a 
polynomial decay of order (r, s) for constants c and h^. 

We note that, if the shape functions are invariant under a linear transformation, then the left 
hand side of (j25p is also invariant under the same transformation change. 

The next proposition establishes a key concentration property of approximation schemes with 
polynomial decay. 
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Proposition 4 (Shape- function concentration). Let {I,W,P} be an approximation scheme. Sup- 
pose that there exists r > such that P has h-density bounded by r. Suppose, in addition, that the 
approximation scheme has polynomial decay of order (r, s) for constants c > and h, with 2s > d. 
Then, for every 6 > there exists a constant ce > such that 

^ \wa{x)\<e, (26) 

XaeP\B(x,cgh) 

everywhere in Q. 

Proof. For every nonnegative integer t > 1, let Ut{x) be the ring of node points Ut{x) = {xa G 
P: (t — l)h < \xa — x\ < th}. Note that P = W^iUtix). By Proposition [21 there exists a constant 
c' that depends on r and d such that, for any t > 1, the number of node points of Ut{x) is at most 
^Ut{x) < c'f^"^. Since the approximation scheme has polynomial decay of order (r, s) with 2s > d, 
for any integer c^ > 1 we have 

oo oo oo 

J2 l^«(^)l ^ E E l"'-(^)l ^ E C't'-'c((t - 1)2 + l)-« < ^ c'f'-'ct-'' 

Xa£P\B(x,Cgh} t=CgXaeUt{x) t=Cg t = Cg 

OO ^ ' 



t=Cg 

Note that the series Ylt^i c'ct~^~^^^~ ' is finite. In particular, there exists a value co < oo, depend- 
ing on d, T, and 6, such that J2t^cg c'ci"^"^^*"'^) <e. D 

For an n-consistent approximation scheme with sufficiently high polynomial decay, the following 
theorem provides a uniform interpolation error bound. 

Theorem 1 (Uniform interpolation error bound). Let {L,W,P} be an approximation scheme. 
Suppose that: 

i) The approximation scheme is n-consistent, n >0. 

a) There exists r > such that P has h-density bounded by r. 

Hi) The approximation scheme has polynomial decay of order (r, s) with 2s > d + m -\- 1, where 
in = min{n,^}. 

Let u S C^"'"^(conv(J7)). Then, there exists a constant C < oo such that 

\D''ui{x) - D''u{x)\ < C ||L>™+^u||^ /i™+i-l°l, (28) 

for every \a\ < min{?7i,r} and x G il. 

Proof. By Proposition [3l 

|D%/(x) - Z)"n(x)| < Y, \R^+i{xa,x)\ \D''wa{x)\ (29) 

for every multiindex a of degree less or equal to min{m, r} and every x ^ Vt. Next, we proceed 
to bound the right-hand side of this inequality. For each nonnegative integer i > 1, let Ut{x) be 



the ring of nodal points Ut{x) = {xa G P: {t — l)h < \xa — x\ < th}. Note that P = W^iUtix). 
By Proposition [21 there exists a constant c that depends on r and d such that, for any t > 1, the 
number of node points of Ut{x) is at most ^Ut{x) < ct'^~^. In addition, from (j20p we have 



\Rm+liXa,x)\ < 



d 



rn+l 



D'^^'uW (th) 



m+1 



(m + 1)! 

By the assumption of polynomial decay there exists a constant < c' < oo such that 

2 \ -^ 



\D''Wa{x)\ < C' 






h 



for every Xa G Ut{x). From the preceding bounds we have 

^ \Rm+l{Xa,x)\ \D°'Wa{x)\ 



ael 



^ ^ \Rm+l{Xa,x)\\D°'Waix)\ 
t=l Xa£Ut{x) 



oo 

<T T 

~ ^-' ^^ (m + l)! 






(m + l ! " "°° ^-^^ ^ ' 



m+l-2s 



t=l 

oo 



fm + D! " "°° ^ 



m+l-2s 



t=l 

oo 



Since d + m - 2s < -1, it follows that YT=\ i^+'^+^-2'* < oo. Thus, 

I -(V/ v/i — no ' 



(30) 



+ 1 /i-l"l <c'((i- 1)2 + 1) '/i-l"l <5c't-2^/i-l"l, (31) 



(32) 



(33) 



for every x G fi, where we note that the constant C = hdc-, — xry (X^^^ ^d+m 2s^ depends on r, d, 



c', and s. 



D 



The following corollaries to Theorem [T] show that a function in W^'^^^t) can be approximated 
by means of consistent approximation schemes of polynomial decay. 



Corollary 1. Under the assumptions of Theorem\^ 



F-f ^lwi.p(n) 



< C7||Z)'"+^n||„o/i^+"'"^ 



(34) 



for 1 <p < oo, j = min{n, r, £} and every u € C '^^(conv(r2)). 



Proof. By Theorem [H there exists a constant < C < oo such that, 

hi - uWc^in) < C\\D"'+\\\^h'+"'-^ , (35) 

so that the assertion follows from the continuous embedding C^{Cl) ^^ W^'P{i}). D 

Convergence in W^'^{i}) finally follows from standard theory of approximation by continuous 
functions (cf. e. g., [T]). For completeness, we proceed to note a particular case of practical relevance. 
We recall that a domain Q satisfies the segment condition if, for all x in the boundary of Q, there 
exists a neighborhood Ux and a direction yx ^ such that, for any point z & Cl n Ux, the point 
z + tyx belongs to Q, for all < t < 1 ([1], §3.21). Convex domains satisfy the segment condition 
without additional restrictions on their boundaries. We additionally recall ([1], thm. 3.22) that, if il 
satisfies the segment condition, then the functions of C^(M ) restricted to Q are dense in W^'^{Q,) 
for 1 < p < oo. 

Corollary 2. Let il. be a domain satisfying the segment condition. Suppose that the assumptions 
of TheoremUl hold uniformly for hk -^ 0. Let 1 < p < oo and j = min{n, r}. Then, for every 
u € W^'P{n) there exists a sequence Uk G Xk such that u^ — )• u. 



Proof. By the density of C^(M ) in W^'P{Q), there is a sequence of functions Vi G C^(M ) whose 
restrictions to il. converge to u in W^'^{il.). The corollary then follows by approximating each Vi by 
a sequence Uj^. S Xk and passing to a diagonal sequence. D 

4 Application to the Local Maximum-Entropy (LME) Approxi- 
mation Scheme: Interior Estimates 

In this Section, we specialize the results of Section [3] to the LME approximation schemes. We begin 
with a brief review of the definition and some of the properties of the Local Max-Ent Approximation 
scheme of Arroyo and Ortiz [2] (see also O [27j for a description of the method, and [24^ [25} 
[26] for related work). We recall that a convex approximation scheme is a first-order consistent 
approximation scheme {/, W, P} whose shape functions are non-negative. Convex approximation 
schemes satisfy a weak Kronecker-delta property at the boundary (cf. [2]), i. e., the approximation 
on the boundary of the domain does not depend on the nodal data over the interior points. This 
property simplifies the enforcement of essential boundary conditions. As pointed out in |2], in a 
convex approximation scheme the shape functions Wa{x), a a I, are well-defined if and only if 
X G conv(i-*). Therefore, for such schemes to be feasible the domain must be a subset of conv(P). 
The Local Maximum-Entropy (LME) approximation scheme [2] is a convex approximation 
scheme that aims to satisfy two objectives simultaneously: 

1. Unbiased statistical inference based on the nodal data. 

2. Shape functions of least width. 

Since for each point x, the shape functions of a convex approximation scheme are nonnegative 
and add up to 1, they can be thought of as the probability distribution of a random variable. 
The statistical inference of the shape functions is then measured by the entropy of the associated 



probability distribution, as defined in information theory |22 1 I14 1 [T5]. The entropy of a probabihty 
distribution p over / is: 

H{p) = -"^Pa^ogpa, (36) 

where log = 0. The least biased probability distribution p is that which maximizes the entropy. 
In addition, the width of a non-negative function w about a point ^ is identified with the second 
moment 

U.{w)= f w{x)\x-i\^dx. (37) 



Jn 



Thus, the width U^{w) measures how concentrated w is about ^. According to this measure of 
width, the most local approximation scheme is that which minimizes the total width 

U{W) = Y,UaiWa)= j Y.Wa{x)\x - Xa\^ dx. (38) 

a&I -^^ aG/ 

The Local Maximum-Entropy approximation schemes combine the functionals (j36p and (|38p into 
a single objective. More precisely, for a parameter /3 > 0, the LME approximation scheme is the 
minimizer of the functional 

Fp{W) = fiU{W) - H{W) (39) 

under the restriction of first-order consistency. Because of the local nature of this functional, it can 
be minimized pointwise, leading to the local convex minimization problem: 



mm.fj3{x,w{x)) = ^^Wa{x)\x - Xa\^ + -y^ Wa{x) log Wa{x)^ 

subject to: Wa{x) > 0, a (^ I, y^Wa{x) = 1, ^ ^ Wa{x)Xa 

a&I a^I 



(LME) 



In the limit of /3 ^- oo the function fp reduces to the power function of Raj an [19] . whose minimizers 
define the piecewise-affine shape functions supported by the Delaunay triangulations associated with 
P. 

Next we collect alternative characterizations of the LME shape functions based on duality 
theory. Let Z : M*^ x M'^ ^^ R be the partition function 

Z{X, A) = ^ ^-P\x-Xa\^ + {X(X),X-Xa) (40) 

ael 

of the point set. For every point x G conv(P), the problem (LME) has a unique solution {w*{x) : a G 
/}. Moreover, for every point x € conv(P), the optimal shape functions w*{x) at x are of the form 

g— ;S|x— Xap+(A*(x),a;— Xa) 

where the vector X*{x) € M'^ minimizes the function 

log Z(X, A) = log ( ^ g-/3|x-x,p + (A,x-x,) j _ (42) 

\ael J 
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At points X belonging to the boundary of convP, the shape functions take expressions similar 
to ()4ip that solely involve the node points on the minimal face of conv(P) that contains x. The 
gradient of logZ{x, A) with respect to A is 

d 

r{x, A) = ^ log Z{x, A) = ^ Waix, A)(x - Xa). (43) 

In addition, the Hessian of log Z{x, A) with respect to A follows as 

J{x, A) E^ --2 log Z{x, A) = ^ Wa{x, \){x - Xa) ® {x - Xa) - r{x, A) (g) r(x, A). (44) 

ael 



Since r{x,X*{x)) = 0, 



J*{x) = J{x,X*{x)) = ^Wl{x){x - Xa) <^ {X - Xa). (45) 



It can be shown that J*{x) is positive definite. In addition, the optimal shape functions w^ : conv(P) 
M are C°° and have gradient 

Vwlix) = -w:ix)iJ*ix))-Hx - Xa). (46) 

We refer the reader to [2] for the proofs of the preceding results and identities. 

The following lemma shows that, for a point set P that is an /i-covering of its closed convex 
hull conv(P), and for every point x G conv(P), for any vector A ^ there exists at least one node 
point xx in P that is close to x, and such that x — xa is closely aligned with A. 

Lemma 2. Let P be a finite point set that is an h-covering of its convex hull conv(P) for some 
h > 0. Let X be a point in conv(P). Let e > be such that B{x,eh) C conv(i-'). Let A 7^ € M"^. 
Then, there exists a node point xx ^ P such that 

i) eh < \x — xx\ < (e + 1)^, 

ii) \X\eh < (A, X — xa). 

Proof. Let x be the point x = x — ^A. Since the distance between x and x is eh, x G B{x,eh) C 
conv(P). In particular, since the point set P is an /i-covering of conv(P), there exists a d-simplex 
Tx of size at most h, with vertices in P that contains x. Let H^ C M'^ be the halfspace {z G 
M : (A,x — z) > 0}. The point x belongs to H~^. Moreover, since the d-simplex Tx contains x, it 
follows that at least one extreme point of Tx also belongs to H~^. Let xa be that extreme point. 
Note that xa is also a node point of P. We have the estimate 

l^; — xa| < |x — x| + |x — xa| < eh + h = {e + l)h. (47) 

In addition, we have 

|x — xaP = |x — xp + |x — xaP + 2(x — x,x — Xa). (48) 

By the definition of x, and since xa belongs to H^ , it follows that 

eh 
{x-x,x- Xx) = j-r{X,x-xx) >0. (49) 
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From this inequality and eq. ()48p we obtain 

|x-XAp > |x-xp = (e/i)2, (50) 

or Ix — xaJ > eh. Finally, from the definition of x we have 

{X,x-x) = \X\£h (51) 

and 

(A, X — xx) = {X, X — x) + (A, X — x\) > \X\eh, (52) 

where we have used that xx belongs to if"*" and, hence, (A, x — x\) > 0. D 



In view of (14ip . in order to verify that the LME shape functions have polynomial decay we 
require a bound on the minimizer X*{x) € M"' of the partition function Z{x, A), eqs. ()40p and ()42p . 
To this end, we begin with the following lemma. 

Lemma 3. Let P be a point set that is an h-covering of il. with h-density bounded by r, for some 
h, T > 0. Let j3 = jj^ for some 7 > 0. Then, there exists a constant cz that depends on 7, t, and 
d, such that 

Z{x,0)<cz (53) 

for every x G conv(P). 

Proof. For every nonnegative integer t > 1, let Ut{x) be the subset of node points Ut{x) = {xa G 
P: {t — l)h < \xa — x\ < th}. Then, by Proposition [2] we have 

00 / 
Z{x, 0) = ^ e-^l^-^"!' = ^ Y. e-^'""""'' 



< f2 (#C/t(x)e-^(*-i)''^') < X;ct"-^e-^(*-i)' = cz. 



t=i t=i 

It is readily verified that the series of the right hand side is absolutely convergent. Moreover, 
because this series is defined in terms of 7, r, and d, its limit cz also depends on 7, r, and d 
only. n 

By optimality, A*(x) has the property that Z{x,X*{x)) < Z(x,0). This observation, combined 
with the upper bound on Z(x,0) of Lemma [31 suffices to estimate |A*(x)|. 

Lemma 4. Let P be a point set that is an h-covering of Q with h-density bounded by t, for some 
h, T > 0. Let /3 = ^ for some 7 > and e > 0. Then, there exists a constant cx > that depends 
on 7, T, and d only such that 

\y{x)\ < . f\., (55) 

for every point x such that B{x,eh) C conv(P). 
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Proof. We note that, since log is an increasing function, X*{x) also minimizes Z{x,X). Let 82 = 
min{e, 1}. We proceed to find a constant cx such that, if |A| > f^, then Z{x,X) > Z{x,0). To 
this end, let A 7^ be a fixed vector. Since B{x,e2h) C conv(P) and since P is an /i-covering 
of conv(P), by Lemma [2] there exists a point x\ ^ P such that £2^ ^ \x — x\\ < (^2 + ^)h and 
{X,x — xx) > |A|e2^- Using these inequalities and noting that £2 < 1, we further obtain 

By Lemma [3l there exists a constant c^ that depends on 7, r, and d, such that Z(x,0) < cz- 
Combining this bound with eq. (|56p. it follows that a sufficient condition for A not to be optimal is 
that e~^'^"'"^2'i|A| y (.^ or, equivalently, 

^ Inez + 47 ^ Cx ,^^. 

min{e, l}/i min{e, l}/i 

Therefore, (155p is a necessary condition for X*{x) to be optimal. D 

We note that, for fixed e > and for points x at distance e or greater to the boundary of 
conv(P), the upper bound ()55p is 0{h~^). By contrast, for points x € conv(P) arbitrarily close of 
the boundary of conv(P), the right hand side of ()55p diverges. The following example shows that 
|A*(x)| may indeed diverge near the boundary. 

Example 1. Let 0, = [a,b] cM., h = b — a and let P = {a, b} be a point set of Q. Let /3 = ^ for 
a some 7 > 0. The optimality condition for X*{x) is 

^^^ = e-^(-a)^+A'(x){x-a)^^ _ ^^ ^ ^.^^,_a-h)^+x'ix)i.~a-h)^^ _ a - h) = 0. (58) 

For this condition we find 

A*(x) = ^osi^ + h-x)^-log{^-a) + ^(2x - 2a - h). (59) 

For a fixed < e < 1, and for points x £ {a + eh, a + h — eh), we indeed have |A*(x)| = 0{h~^). 
However, lim3,_j.„+ A*(x) = oo, and lim2._s.5- ^*{^) — — cjo. We note that the LME shape functions 
for this case reduce to 

*i \ o- + h-x 
Wa{x) = , (60a) 

<{x) = ^. (60b) 

In particular, the shape functions and their derivatives are bounded in Vt even though the value of 
|A*(x)| is unbounded at the boundary. From a computational perspective, this example suggests 
that computing the shape functions and their derivatives using Equations (P2|) and ([^T]) may be 
unstable near the boundary, even if the shape functions and their derivatives are themselves well- 
behaved. In Section [5] we will examine the behavior of the shape functions near the boundary more 
thoroughly. D 

The following lemma supplies the requisite estimate of the partition function Z. 
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Lemma 5. Under the assumptions of Lemma^ there exist constants inz, Mz > that depend on 
7, T, £, and d only and such that 

mz <Z{x,X*ix)) <Mz (61) 

for every point x such that B{x,£h) C conv(i-*). 

Proof. By optimality, Z{x,X*{x)) < Z{x,0) and, by Lemma [3l Z{x,X*{x)) < cz = Mz for every 
X € conv(P). Since P is an /i-covering of conv(P), there exists a point xq S -P at distance to x less 
or equal to h. In addition, by Lemma HI there exists a constant cx such that |A*(2;)| < -^, where 
£2 = min{e, 1}. We thus have 

Z{x, A*(x)) > e-/5|^-^ol'+<^*W'^-^o> > g-7-|A* Wl k-^ol > g-7-|| = mz > 0, (62) 

as advertised. D 

Recah that J*{x) G M'^^'^ is the Hessian of logZ(x, A*(x)) with respect to A, eq. ([15]). We 
proceed to estimate || J*(x)~-'^||. 

Lemma 6. Let P be a point set that is an h-covering of Q with h-density bounded by r, for some 
h, T > 0. Let /3 = p- for some 7 > 0. Let e > 0. Let x be such that B{x,£h) C conv(P). Then, 
there exists a constant cj-i > that depends on r, 7, £, and d such that 

II j*t \-l\\ — r^ V-^) \y)\ ^ 1,-2 iao\ 

||J (x) II = sup — <Cj-\h . (63) 

2/^0 \y\ 

Proof. Let u 7^ be a fixed vector. Then, from eq. (|45|) we have 

V^ „-l3\x-Xa\'^ + {\*{x),X-Xa) /„, ^_^ \2 

Next, we analyze the numerator and denominator of the right-hand side in turn. Let £2 = min{e, 1}. 
By LemmaEl there exists a point Xu € P such that £2^1 < l^;— x^l < {e2+l)h and {u,x—Xu) > \u\e2h. 
Since B{x,£h) C conv(P), by LemmalHthere exists a constant c\ such that |A*(x)| < ^, and we 
have 

|(A*(x),x - x,)| < |A*(x)| |x - x„| < ^(£2 + l)h = ^(£2 + 1). (65) 

£2/1 £2 

Hence, 

y^ g-/3|2:-a:aP + (A*(a:),a;-a::a>/^ 2; _ 2; \2 

aei (66) 

where we write /3 = ^. Combining the bound supplied by Lemma [5] with eq. ()66p . we get 

2 

I T T*/ \ I \ -(e2+l)^7-7T-(s2+l) £2 I |2t,2 i |272 /c'7^ 

|u J (x)m| > e ' ' £2 ' '^^l^^l^/i^ = cj|ii| n , (67) 

where cj = e '"a ^^^ -^-^ > depends on 7, r, £, and d only. Let Amin(3;) be the smallest 

eigenvalue of J*(x). Since J*{x) is positive-definite [2], it follows that Amin(2;) > 0. Inequality 



(j67|) then implies that Amin(2;) > cj/i^. Since ||J*(x) ^|| = 1/Amm(a^), the estimate (j63|) follows 
immediately with Cj-i = 1/cj. D 



14 



We are finally in a position to estimate the derivatives of the LME shape functions. 

Proposition 5. Let P be a point set that is an h-covering of ft with h-density bounded by r, for 
some h, T > 0. Let f3 = -r^, for some 7 > and e > 0. Let W = {u;* : a £ 1} be the optimal shape 
functions of the LME approximation scheme with node set P and parameter (3. Then, 

|V'u;*(x)| < cj-i wl{x)\x - Xa\ h~^, (68) 

for every point x such that B{x, eh) C conv(i-') and every point Xa € P- 

Proof. The estimate (j68p follows immediately from Lemma [6] and eq. (j46p . D 

Next we show that the LME approximation scheme has polynomial decay of order (l,s) for 
every s > 1. 

Proposition 6. Let {I, W, P} be an LME approximation scheme. Suppose that P is an h-covering 
of Q, P has h-density bounded by t, j3 = '^jh? for some 7 > 0. Let e > 0, and s > 1. Then, there 
exists a constant c > (depending on d, 7, t, e, and s) such that the approximation scheme has 
polynomial decay of order (1, s) for c and h in Q^/,, = {x G M"^ s. t. B{x, eh) C conv(i-')}. 

Proof. We recall that the LME shape functions are C°° on conv(P) ([2]). Next, we show that there 
exists a constant c > that depends on 7, r, d, e, and s, such that, for any A;, 



sup sup sup 1 + 






/il"l|Z)°u;:(x)| <C. 



(69) 



FroixL Lemmas m and [5] we have 



< <(x) 



'■\x-Xa\^ + {X*{x),X — Xa) 

Z{x,\*ix)) 



^~-y\{x-Xa)/h\^+\X*{x)\\x-Xa\ ^-■y\{x-Xa)/h\^+cx\{x-Xa)/h\ 

< < , 



(70) 



mz 



with cx = cx/ i'nm{e, 1}. In addition, 



sup sup 1 + 



< sup sup I 1 + 

xGfl^h a.Gl 



<c :=sup(l+t2)' 
t>0 



X 


•^a 


2\ 


h 


/ 


+ 


X X(i 


h 





mz 

<(x) 
'\ g-'y\{x-Xa)/h\'^+cx\{x-Xa)/h\ 



mz 



-^It^+cxt 



< 00, 



mz 



(71) 



since 



mz 



is a rapidly decreasing function of t. We note that c' is defined in terms of the 



15 



constants d, 7, r, e, and s. Thus, by Proposition [5] we have 

2^ 



sup sup 



< sup sup 



< sup sup I 1 + 




h\Vwl{x) 



hcj-iw^{x)\x — Xa\h 
cj-iwl{x) < Cj-ic', 



as advertised. 



(72) 



D 



The next Theorem bounds uniformly the error of the approximate function Uk and its derivatives 
to a smooth function u and its derivatives. The result is based on Theorem [1] that holds for a general 
approximation scheme. 

Theorem 2. Under the assumptions of Proposition\^ let u € C^(il). Then, there exists a constant 
C > 0, that depends on 7, r, e, and d only, such that 



|D°U7(x) - D''u{x)\ < C||Z?2u||oo/l^"'"l, 

for X e ^eh, |a| < 1- 

Proof. The theorem follows from Proposition [6] and Theorem [TJ 



(73) 



D 



Finally, we are in a position to show that LME approximation spaces on a domain Q.' are dense 
in W^'^{Q) for subdomains C 0' which are compactly contained in U! . This result is derived from 
the polynomial decay of LME schemes, and the density of approximation schemes of Corollary [2l 

Corollary 3. Let Q be a domain satisfying the segment condition, and let Q' be an auxiliary 
domain such that C 0'. Let {Lk,Wk,Pk} be a sequence of LME approximation schemes in il.' . 
Suppose that the assumptions of Proposition \^ hold for {Ik,Wk,Pk} in ^' uniformly for hk — > 0. 
Let 1 < p < 00. Then, for every u € W^''P{0,) there exists a sequence Uk € Xk such that Uk\n — >■ u. 

Proof As n C Q', there exists r > such that, \Jx^nB{x,r) = {x G M'': dist(x,il) < r} C ^'. The 
sequence of approximation schemes {Ik,Wk,Pk}, when restricted to il., has uniform polynomial 
decay (l,s) for any fixed s. Then, the theorem follows from Corollary [2j D 

Corollary [3] guarantees the density of the LME approximates on W^'^{^), provided that the se- 
quence of LME approximation schemes {Ik, Wk, Pk} is defined on a bigger domain fi'. We note that, 
in this case, the LME scheme does not obey the weak Kronecker-delta property at the boundary of 
Q, making it less straightforward to enforce boundary conditions on 17. However, imposing bound- 
ary conditions can be done in this case by using standard Lagrangian multipliers, see e. g. [131 [8]. 
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5 Application to the Local Maximum-Entropy (LME) Approxi- 
mation Scheme: Estimates up to the Boundary 

In Section|4]we have seen that, for a sequence {Ik, Wk, Pk} of LME approximation schemes, we have 
density of the approximation space Xk in W^^^{il). In order to treat boundary value problems, 
however, we need density results up to the boundary of Q. A way to guarantee the density in 
Vl^^'P(r2) is to work with a sequence {Ik,Wk,Pk} defined on a (strictly) bigger domain O', as 
discussed in Corollary [3l In this section, we analyze the density of the approximation space Xk 
when the domain of the LME scheme is 0,. 

While we will see that density can be extended to Wq '^{^) in general, a major technical difficulty 
with estimates up to the boundary comes from the fact that X*{x) blows up as x approaches di^. 
This blowup is indeed a manifestation of the weak Kronecker-delta property at the boundary, as 
A* will blow up in such a way that in the limit no weight is given to nodal data in the interior of fi. 
For general i}, this behavior can become very complicated and lead to blow up of the gradients of 
the optimal shape functions Vwa, with the result that the general convergence scheme of Section [3] 
is no longer applicable. Therefore, for simplicity we restrict attention to the class of polyhedral 
domains. Under generic assumptions, we shall obtain sufficiently strong estimates on Vtt;* near 
flat pieces of the boundary dQ permitting to show that, away from a small singular part of the 
boundary, Sobolev functions can be approximated by linear combinations of shape functions in the 
limit of /i — )• 0. The singular boundary is of finite 2-capacity. With the help of a capacity argument 
we can then establish approximation results with truncated LME functions in H^{Q) for spatial 
dimension d > 2. 

More precisely, in this section we will assume that O is a convex polytope in M'*, P is an h- 
covering for with convP = 0, such that there exists a constant rj > such that {x £ P: < 
dist(a;, 90) < r]h} = 0. Note that then P n dO, is an /i-covering for d^. 

Assume that A = H D dO,, H some hyperplane, is a flat {d — l)-dimensional subset of the 
boundary of f2. With the aim to control Vw*{x) for x in the vicinity of A, our first task will be 
to exactly estimate the behavior of J*{x) in this regime. First note that with a proper choice of 
the coordinate system we may assume that H = {xi = 0} = {0} x R with Q. n {xi > 0} = 0. 
Accordingly, we write 

\*{x) = {X\{x),\'{x)) G M X R"^-! (74) 

and, for x = (xi,x'), 

Z = y^ g-/3|x-x,|2 + (A',x'-0+(x-x„)i A* ^ ^-^5^) 

as/ 

We fix (5 > and consider points x £ Q. with x = {p, x') G M x M. '^ for p small such that 
-65/1(0, x')r]H = Bsfi{0, x') n A. In the following lemmas we will also set h = 1 for arbitrarily large 
0, and recover the general case in Proposition [7] by rescaling afterwards. Generic positive constants, 
denoted c, c', c" or C, C", will be independent of p and the size of O. 

Lemma 7. There is a constant C > such that 

\X'{x)\ < C and Xl{x) > -C. (76) 

Proof. This result follows along the same lines as the proof of Lemma [4] for the boundedness of A* 
in the interior of 0. D 
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In order to investigate Z we split the sum as 

Q.G/A aGl\lA 

where I a collects those indices a for which Xa & A. 

Lemma 8. As p tends to 0, A^ ^- oo such that p\\ — )• 0. 

Proof. Writing Z as in (|77p and noting that {xa — x)i > i] for a ^ Ia and /) > 0, we see that 

oG/a 

and that this lower bound is in fact achieved only if p\\ — )• and \\ -^ oo. D 

In particular, we see that Z still remains bounded from above and from below by positive 
constants. 

In order to estimate J*, we first observe that the optimality condition ^^ = implies 

p y^ g-/3|a:-a:aP(A',a;'-0+pA* ^ V^ ^-l3\x-Xa\'^ + {\' ,x' -x'J-ixa-x)^^ /^^ _ ^\_ ^-j^gN 

a&lA aeI\lA 

Lemma 9. There is a constant c > such that the first entry J^-^ in J*{x) satisfies 

Jii > cp. (80) 

Proof. Since (xq — x)i > r] for a ^ Ia we find by eq. (j79p and Lemma [7] 

> ^-1 y^ g-/3|x-Zap + (A',x'-X^)-(Xa-x)lA*/-^ — x)^ 

a<aI\lA 

>riZ-' Y, e"^l"-"'"l'+<^''"'-<>-("'^-")i^i(xa-x)i (81) 

aG/\/A 

= rjpZ^^ y e-/3|=^-^aP+(A',x'-x'J+pA* 

ael A 

> cp 

for some constant c > 0. D 

We now derive an upper bound for the entries of the first row and column of J*. 
Lemma 10. For any < fi < 1 there exists a constant C > such that 

\Jij\ = \J*i\<Cp^', j = l,...,d. (82) 
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Proof. For j = 1, . . . ,d we have Jj*- = J*i = J2aei ^a(^)(^ ~ Xa)i{x — Xa)j- First summing over 
a € /a gives the obvious bound 



y^ wl{x){x - Xa)l{x - Xa)j 



< 



2-'E 



X~Xar + {\' ,x' -x'^)+p\l 



P\X - Xa\ 



aelA 



(83) 



<Cp. 



In order to estimate the remaining sum, we let p = - and choose 1 < g < oo with - + - = 1, 
so that 



< I Yl w:ix)ixa-xrA Yl <(^)i 

.a€l\lA } \a&I\lA 



y^ Wl{x){x -Xa)l{x-Xa)j 
aeI\lA 

by Holder's inequality. Here the second factor in (j8^ is bounded by 

-/3\x-Xa\'^ + {X' ,x'-x'J-{Xa-x)lXl 






(84) 






3^ ~ 3;a|„- I ^ C, 



(85) 



see Lemma [71 To estimate the first factor we note that, since P is a 1-covering of 0,, there exists 
a & I \ Ia such that 

\X(i X *^ o 

for a constant C' > ry. For p sufficiently small and thus A^ sufficiently large, we then have the 
estimate 

Y wl{x){xa - x){ < Z~^ Y. e-^l"-"'^l'+<^''"'-<>-^'^i(x,-x)? 



aeI\lA 

{Xa-x)l>C' 



aeI\lA 

{Xa-x)l>C' 



(86) 



<CZ-^i]-^ Y e-'^l^-^'"l'+<^*'^-^''>(xa-x)i, 

a&I\lA 

as (xq — x)i > ?7 for a € / \ Ia- On the other hand, for a with (xa — x)i < C we have the bound 

{Xa-x)l<C{Xa-x)i. (87) 

Combining the two last estimates, we see that the term in the first factor of (184p satisfies 

Y '^*ai^)i^ - ^afl < C Y Wl{x){x - Xa)l. 
a&I\lA aeI\lA 

Since by (j79p this last expression is bounded by Cp, we arrive at 



Y wl{x){x - Xa)l{x - Xa)j 



aeI\lA 



< Cpp = Cp'' 



(89) 
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by (j84p . Together with the bound (j83p for the first part of the sum we have shown that indeed 

\Ju\<Cif. 

D 

For the remaining part B = {J*j)2<i,j<n of the matrix J* we obtain the following lower matrix 
bound. 



Lemma 11. There is a constant c > such that 

B>cldn-1. 



(90) 



Proof. As P n if is a 1-covering for A, there is a set J = {oi, . . . , a^-i} C /^ of d — 1 points such 
that c' < \x' — x'^\ < c" and det(x' — rr^^ , . . . ,x' — x'g^^_-^ ) > c' for suitable constants c' and c" . Then 



> z-1 ^e-^l^-^'^l'+<^''^''-^^>+''^i(x' - x^) (» (x' - x'J 

aeJ 

> c^{x' - x'J (g) {x' - x'^) 



(91) 



> cldd-1 
since all the projections {x' — x^) ® (x' — x^) are nonnegative. D 

As a consequence of the above results, we obtain an estimate for the inverse matrix (J*)~^ = 

Lemma 12. For any < fi < -^ there exists a constant C such that 

'<Cp-^ fori=j = l, 

<Cp-^' fori = l,j = 2,...,dorj = l,i = 2,...,d and (92) 

<C fori,j = 2,...,d. 



\Jij\ \ 



Proof. First note that, expanding with respect to the first row, for ^ < fi < 1 we have 

det J* = Ji*i det B + 0{p^^') = J^ det B > cJ^^ > cp 
by Lemmas m [To] and [TTJ Furthermore, as \J*\ < C, we have 



(93) 



|(cof J* 



Hj\ 



= |det5|<C fori = j = l, 

< Cp^ for i = 1, J > 2 or J = 1, i > 2 and 

[<C{J*n+p^n fori,i>2. 



(94) 



for C sufficiently large. Now, Cramer's rule 

(J*)-i = (detJ*)"i(cof J*)^ 



(95) 
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implies 



\Jij\< 



< Cp ^ for i = j = 1, 

< Cp^-i for i = 1, j > 2 or j = 1, i > 2 and (96) 
<C for i,j> 2 



and thus the assertion follows by choosing /x such that fi = 1 — jl. D 

Lemma 13. For any s > and < /i < 2 there is a constant C > such that 

{1 + \X- Xa\y |V<(X)| < C{1 + p-f'lx - Xa\). (97) 

Proof. If a € /a ) then x — Xa = {p, x' — x'^) and Lemma [T2] shows 

I {rr\x -Xa)\< C(l + p-I'lx' - <|) < C{1 + p-I'lx -Xa\). (98) 

So, by dlBD, 

|V<(x)| < C\w:{x)\{l + p-'^lx - Xa\). (99) 

Now, using that (1 + |a; — a;ap)''|tt'a(^)l — ^ fo^ ^^Y '^j "^^ see that the estimate holds true for 
a G Ia- 

On the other hand, if a ^ J^i, then Lemma [T2] only gives 

\ir)-\x-Xa)\<Cp'^\x-Xa\, (100) 

whence 

|V<(x)| < C7K(x)|p-V - Xa\. (101) 

But since (x^ — x)i > rj for all b ^ Ia, ^^ also get 

(l+|x-xj2)«|-i^*(2;)| 






(102) 



66/\/a 



This term can now be estimated by Cp^p ^ for < /i < 1 precisely as the left hand side of ()84|) in 
Lemma [TOl which leads to 



(1 + |X - Xa\'^y\Vwl{x)\ < Cp-^'lx - Xa\ (103) 

for < /i < 1. D 

Undoing the rescaling of h we can now summarize the previous lemmas in the following propo- 
sition the boundary behavior of Vw* near flat parts of 90. 
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Proposition 7. Suppose x = {p, x') gM. xM.'^ ^ is such that Bsh{0, x') D H = Bsh{0, x') n A for a 
boundary (d — l)-face A = d^ n H. Let s > and < /i < 2 • There is a constant C > such that 



1 + 



X Jvn 



h 



h\Vwl^,^{x)\ <C{l + hf'dist-''{x,dn)) 



(104) 



Proof. If P is an /i-covering for il, then h ^P is a 1-covering for h ^Q. Using subscripts h to 
highlight the dependence on h, we have 



Zh{x) = J2 ^ 



-^\x-Xh\'^ + {X*^^,x-xi,) 



Xb€P 



y2 g-^i^i'+^'^^w''^^ 



Xb&P 



(105) 



y^ e-7lf-^bl'+('^A*^).f-^b>_ 



This expression is minimized at hXIf^Jx) = A/j^n(|). For the shape functions t/;* ^ we denote by Wa,i 
the shape function corresponding to the node ^ € h~^P and obtain 
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Applying Lemma [13] to uj* i we therefore obtain 

2\ " / 
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From our previous calculations we know already that the left hand side is bounded by a constant 
away from the boundary of 0,. Suppose x is such that dist~'^(^,^) > 1. Since s is arbitrary we 
can absorb the last factor on the right hand side into the prefactor of the left hand side and finally 
obtain 



1 + 






h\Vwlh{x)\ <C(l + /i'^dist-^(x,5fi)) 



(108) 
D 
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Now suppose that x is a general point near a possibly lower dimensional edge of dil. More 
precisely, x is close to an ?7i-face of A of dQ, which is the intersection oi d — m hyperplanes 
Hi, . . . , Hd-m with linearly independent normals which constitute dVt in the vicinity of x. 

Lemma 14. There exists R > such that for all Xa & P with dist(2;(j, dQ) > Rh 

2\ * 



1 + 






h 



h\Vw:Ax)\ < 1. (109) 



Proof. We first assume again that h = 1. Let H be the hyperplane containing x which is perpen- 
dicular to A*. Similarly as in Lemma [7] we see that, as dist(x,^) -^ 0, |A*| tends to infinity such 
that the projection of A* onto P| Hi remains bounded and that there are constants c, C > such 
that 

y-x,—\ =dist{y-x,H) >cdist{y,HiU...UHd-m)-C (110) 

for every y £ Q. 

Choose a set J = {ai, . . . , a^} C / of d points such that c' < \x — Xa\ < c" and det(x — 
Xai , . . . ,x — Xa^) > c' for Suitable constants c' and c". Then 

J* >Z-^^e-^l^-^'^l'+<^*'^-^">(^-2;a)®(x-Xa) >ce-^'^*lldrf_i (111) 

since all the projections {x — Xg) ® {x — Xa) are nonnegative. It follows that 

(J*)~i <Ce^l^*lldd_i. (112) 

Now if Xa G il satisfies dist(2;, 5il) > R, then 

{x-xa,\*) < -{cd[st{xa,HiU...UHd-m)+C)\X*\ < {-cR + C)\X*\. (113) 

So 

\w*\ < ^■"^e^^l^"'^"!^^^'*'*'''"'^''^ < ^-ig"'^!'^"^''!^"^''-'^"'^^!^*'. (114) 

It follows from (1461) that 



l + \x-xa\') \Vwl{x)\ < C (l + \x - xa\^j g"/?!^-^.! -(^«+^)l^*le^l^*l |x - x„ 

< Ce^^C-cR)\X'\ < ^ 



(115) 



for R sufficiently large, which proves the Lemma for h = 1. The estimate for general h now follows 
directly by rescaling as before. D 

We are now in a position to prove our main density results up to the boundary. Density in 
Wq'^{0,) in fact only relies on our previous interior estimates, see Section [U and Lemma HM and is 
true for general, not necessarily polyhedral domains O. 

Theorem 3. Let 0, be a bounded polyhedron and {Ik,Wk,Pk} be a sequence of LME approximation 
schemes satisfying the assumptions of Proposition uniformly for h^ — t- 0. Then for any u G 
Wq'^{VI), 1 <p < oo, there exists a sequence Uj, G Xj, such that Uk — > u. 
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Proof. It suffices to consider u G C^(0). Let Uk = uj^,. By Proposition [3] we have 

\D''uk{x) - D''u{x)\ < Y, \R2{xa,x)\D''wl{x)l (116) 

which for dist(x, 5il) > eh can be estimated by C||D^u||oo^^. — >■ as /i^ — >■ 0, see Theorem[TJ 

If dist(x,5ri) < eh, then with the help of Lemma [TJ] precisely the same arguments as in the 
proof of Theorem [1] show that 

Y, \R2{xa,x)\D^wl{x)\<C\\DMoohl~^''^^Q (117) 

a&I 

\x — Xa\>Rh 

for a sufficiently large constant i? > 0. But the remaining part of the sum vanishes for small 
/ifc as R2{xa,x) = if Ix — Xa\ < Rh, since then u vanishes on a neighborhood of the segment 
{xa + e{x -Xa): 9 e [0, 1]} C {x € Q : dist(x, dn) <{R + e)hk}. D 

In order to formulate our main result on density up to the boundary we denote by d^O, the 
union of the interiors of the {d — l)-faces of dQ. (5*0 is the reduced boundary in the language of 
geometric measure theory.) The part of O a distance eh away from the singular boundary 90 \d^:Q 
is denoted Og/j = {x £ Q: dist(x,5r2 \ 9*0) > eh}. 

Theorem 4. Let 0, be a bounded polyhedron and {Ik,Wk,Pk} be a sequence of LME approximation 
schemes satisfying the assumptions of Proposition^^ uniformly for h^ — s- 0. Let e > Q. Then for 
any u G W^''P{Q), I < p < oo, there exists a sequence Uk G X^ such that \\uk — u\\^ri..pfQ ) — ^ 0. 

Therefore, density holds away from the singular boundary. 

Proof. Let u G C^(0). As in the proof of Theorem [3] we find by Proposition [3] and the arguments 
in the proof of Theorem [T] that for all x G ^eh 

|D°n/(x)-D"n(x)| < CIlD^nlloo/i^"'"' {l + h^'dist->'{x,^n)), (118) 

< ^ < i, where in addition we have applied Proposition [7] in order to estimate Vw* near the 
regular boundary. Consequently, 

/ |L>"u/(x) - D'^u{x)\Pdx < C/i(2-l"l)p / 1 + hPf'dist-Pf' {x,dn) dx 

< C/i(2~l"l)p fi + /jP/^ f t"PM (ii\ (119) 

= C/i(2-|"l)P 

for /i sufficiently close to 0. D 

Note that since the {d — 2)-dimensional Hausdorff measure of 90 \ 9*0 is finite, this set has 
zero 2-capacity for d > 3. Theorem [H thus shows that u can be approximated by Uk G Xk in H^ 
up to sets of arbitrarily small 2-capacity. With the help of a capacity argument we obtain from 
Theorem HI 
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Corollary 4. Let Q be a bounded polyhedron and {Ik, Wk,Pk} be a sequence of LME approximation 
schemes, e > 0. Suppose d > 2. Then for any u S H^{Q) there exists a sequence Xk £ H^ with 
Xk ^ ^ in H^ o,nd a sequence Uk S Xk such that XkUk -^ u in H^. 

Proof. Since the {d—2) dimensional Hausdorff measure of the singular part dft\dt:il. of the boundary 
is finite, this set has zero 2-capacity: 

Cap2{dn\dM) = 0. (120) 

In particular, for every neighborhood V of dQ \ 8^,0, and 6 > there exists a function ips S H^{W^) 
with compact support in V such that ^^ > 1 in a smaller neighborhood of dQ \ 9*i7 and 

\ViJs\^dx<6. (121) 

(This follows, e. g., from Theorem 3 and its proof in jlH pp. 155-157].) By replacing, if necessary, 
ipS with a mollification of m.ax{ioam{iJ;s, 1}, —1} we may assume that ips is smooth, \ips\ < 1 and in 
particular tpg = 1 near dQ \d<t:Q. 

Now suppose u G 0^(0,). By Theorem H] we find Uk G Xk with \\uk — 'u||^yi,p/Q ^ ^- for all 
p < oo. Since 1 — ips vanishes in a neighborhood of dQ \ 5*^2, it follows that 



2 

k' 

(122) 



11(1 - 'll^s){Uk - «)llHi(f^) ^ 111 - V'5|lLoo(f2)||ufc - «llj:^i(Q^j + l|V^5||L2(n)lhfc - ^IIloo(Q^J 

by Sobolev embedding. As 



11(1 - ^ps)u - u\\jji = ll^^nll^i 

< / \u\^dx+ / \S/u\'^ dx + \\u\\l^ / iVipsl'^dx (123) 

Jv Jv Jv 

<Ci\V\ + 6) 

and V and 5 can be chosen arbitrarily small, by choosing diagonal sequences we see that every 
u G C^^fl) and hence in fact every u G H^{Cl) can be approximated by sequences (1 — ips^)uk, 

Uk G Xk- n 

6 Concluding remarks 

The preceding analysis shows that, whereas the density of the LME approximating scheme in the 
interior of the domain follows directly from the general results for meshfree approximation schemes, 
the density of the scheme up to the boundary is a matter of considerable delicacy. This situation 
strongly suggests relaxing the positivity constraint and allowing for signed basis functions. This 
relaxation is also required for the formulation of higher-order approximation schemes, as noted by 
[21 [9]. Indeed, in the finite-element limit shape functions of quadratic order and higher are signed 
functions in general. As an additional bonus, signed shape functions enable the consideration of 
general — not necessarily convex — domains. These extensions are pursued in a follow-up publication 
[3], where LME- type approximation schemes of arbitrary order and smoothness are derived and 
their convergence properties are analyzed using the general analysis framework developed in this 
paper. 
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